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' Abstract. In this third paper of our series on near-IR and optical photometry of globular cluster systems in 

(«~^ I early-type galaxies we concentrate on the photometric results for NGC 5846 and NGC 7192, two group ellipticals, 

— , ■ and on a first comparison between the globular cluster systems investigated so far. In NGC 5846 the colour-colour 

' diagram shows clear bi-modality in (V — K), which is confirmed by a KMM test. The mean colour of both peaks 

O-i] were estimated to be {V - A')biue=2.57±0.06 and {V - iC)rod=3.18±0.06. The situation in NGC 7192 is different, 

^ ■ in that the colour-colour diagram gives no evidence for a distinct second population of globular clusters. Using 

^ ' simulated colour distributions of globular cluster systems, we make a first step in quantifying the cumulative age 

, distribution in globular cluster systems. Also here the result for NGC 5846 leads us to the conclusion that its 

' metal-rich globular cluster population contains two globular cluster populations which differ in age by several 

' Gyr. The age structure for NGC 7192 shows instead strong similarity with a single-age population. 
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1. Introduction ity and their age, the discussion about the origin of those 

populations is still ongoing. Starting with Peebles & Dicke 
As shown by many studies during the last decade (^ggg) ^j^^ assumed globular clusters to be the first ob- 
(Zepf & Ashman 1993, Ashman & Zepf 1998, Kundu & -^^^^ ^^^^^^ ^j^^ ^^^^^ universe, the number of possi- 
Whitmore 2001a, b, Larsen et al. 2001, Kissler-Patig et formation scenarios increased drastically since then, 
al. 2002) globular cluster systems are a very powerful ^^^^^ -gg^^g ^^^^^^^ multi-modality 
tool in galaxy formation and evolution studies. Although ^^^^^ colour-colour diagrams of globular cluster sys- 
there is a wide agreement about the existence of sub- ^^^^^ ^^^^^^^ it ^^^^^^ ^^^^ different globular clus- 
populations m cluster systems regardmg their metallic- populations are produced during strong star formation 
events, but the nature of these events remains under de- 
Send offprint requests to: Maren Hempel ^^^^ Besides the merger scenario, favoured by Ashman 
* Based on observations at the Very Large Telescope of the ^ ^epf (1993) (see also |Ashman fc Zepf 1992| Whitmore 
European Southern Observatory, Chile (Program 63.N-0287 Whitmore & Schweizer 1995; Schweizer et 
Based on observations made with the NASA/ESA Hubble , t^. i t-> , • i • i r i 
o rr. , • J A. J i v.- 4- 4-1, c 8.1. 1996; Kisslcr-Fatig 2000), there is a number of al- 
Space telescope, obtained from the data archive at the Space ' o /i 

Telescope Science Institute. STScI is operated by the associa- ternatives, i.e. the accretion scen ario ^(^ote et al. 192 ^1 
tion of Universities for Research in Astronomy, Inc. under the liM£et al. 2002" 'Hilker et al. 1999) or the monolithic col- 
NASA contract NAS 5-26555. lapse (,Forbcs et al. 1997a, .Kisslcr-Patig ct al. 1998| and 
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references therein). From the photometric point of view 
it has been shown that the multi-modahty of the colour 
distribution is a common feature of globular cluster sys- 
tems ( jGebhardt Hi Kissler-PatigTQQQ^ . In particular, the 
existence of a blue, old and metal-poor cluster population 
IjAshman fc Bird IQQSllHurgarella et al. 2001| ) is a general 
feature. Regarding the red sub-population the almost only 
consensus which has been reached so far, is about the large 
varieties between different cluster systems. This includes 
the possible existence of sub-populations within the red 
population which is interpreted as a sign of different star 
formation events in later stages of the galaxy evolution. 

The main problem in specifying and dating dif- 
ferent star formation events arises from the age- 
metallicity degeneracy of the most commonly used 
optical colours. So far the bulk of high-quality pho- 
tometric investigations have been performed using 
Hubble Space Telescope (HST) in the optical wavelength 
regime (|E^hes et a 1 . 1 9981 IKundu fc Whitmnre 1 99SI 
IGebhardt fc Kissler-Patig 19991 iLarsen et al. lO?)!)!. It 
has been shown ( Minniti et al. 19"96l |Kissler-Patig 2000| 
Kissler-P atig et al. 2002, (hereafter cited as Paper I); 
Puzia et al. 2002J (hereafter Paper II)) that combining 
optical and near-infrared data is a more promising 
method to separate age and metallicity effects and to 
access the relative ages of the globular cluster popula- 
tions. The method relies on a sampling effect, where the 
V-band is dominated by stars near the turn-off (TO) 
region, whereas the main contribution to the K-band is 
from giant branch stars l|Yi et al. 2001(1 . Whereas the TO 
is dominated by age effects, the giant branch shows a 
high sensitivity to metallicity IjSaviane et al. 2000|l . This 
results in a similar dependence of {V — I) and {V — K) on 
the age of the clusters, but a higher sensitivity of {V — K) 
to the metallicity. 

In the previous two papers of this series (Papers I and 
II), a systematic survey of globular cluster systems of E 
and SO galaxies in the combined optical and near-infrared 
wavelength range (using the V , I, and Kg bands) has been 
started. In Paper I, we compared the globular cluster sys- 
tems of two ellipticals in the Virgo cluster, namely the 
giant central galaxy M87 and an intermediate-luminosity 
galaxy NGC 4478. We found that in those cases, the 
(V—K) colour distribution yielded roughly consistent con- 
clusions relative to those derived from the {V — I) colours 
measured by HST. In Paper II however, the {V—K) colours 
of globular clusters in NGC 4365 led us to postulate the 
existence of a significant population of intermediate-age 
(^2-6 Gyr old), metal-rich globular clusters, which was 
not revealed by the (V—I) colours. This important result 
has recently been confirmed by deep spectroscopy (Larsen 
et al. 2003), adding credibility to the results derived from 
our optical 4- near-infrared imaging program. 

One of the aims of this series is to study the the glob- 
ular cluster systems of galaxies in different environments, 
e.g. giant ellipticals in centres of clusters as well as rather 
isolated and less luminous galaxies. The present work will 
focus on NGC 5846, a giant EO galaxy in the centre of the 



Table 1. General information about the host galax- 
ies N£C_5846_jmd_J^^^ The references are 
(l) ^de Va ucouleurs et al. l99l, (2): [Schlegel et al. 1998 
(3) Buta fc WiUiams 1995> (4): |Frogel et al. 1978 
(5) T onry et al. 200T1 



Property 


NGC5846 


NGC 7192 


Reference 


RA(J2000) 


15h 06m 29s 


22h 06ni 50s 


(1) 


DEC(J2000) 


+01° 36' 25" 


-64° 18' 57" 


(1) 


Bt,o 


10.87 


12.19 


(1) 


Hb-v 


0.055 


0.034 


(2) 




0.96 


0.92 


(1) 


{V - /)off,o 


1.28±0.01 


1.24±0.01 


(3) 


{V - K),s,o 


3.51±0.01 




(4) 


(m - M)v 


31.98±0.20 


32.89±0.32 


(5) 




-22.07 ± 0.20 


-21.62 ±0.35 


(1),(5) 



Virgo-Libra Cloud, and on NGC 7192, an isolated ellipti- 
cal with only one companion. Basic informations on both 
galaxies are provided in Table 1. 

Both systems have already been studied in the opti- 
cal HForbes et al. 1997bl[Gebhardt fc Kissler-Patig T999| ), 
and show a very broad {V — I) colour distribution. This 
work aims at the detection of different globular cluster 
sub-populations, the determination of their age structure 
and a first comparison of various globular cluster systems. 
The present paper is organised as follows. In §2 the obser- 
vations and the data reduction procedures are described. 
Chapter 3 contains the main results of the observations 
and §4 describes our approach towards quantifying the 
age structure in globular cluster systems and the results 
for both systems. In chapter 5 we will give an outlook on 
the upcoming work. 

2. Observations and Data Reduction 

2.1. VLT/ISAAC Near-Infrared Data 

The Kg band exposures for both galaxies have been taken 
in service mode (ESO program 63.N-0287) with the Near- 
Infrared Spectrometer And Array Camera (ISAAC) at- 
tached to the Unit Telescope 1 (Antu) of the European 
Southern Observatory's Very Large Telescope (VLT). The 
field-of-view of ISAAC'S Rockwell infrared array is 2.5' x 
2.5', with a pixel scale of 0.147"/pixel. All data were ob- 
tained in April and June 1999 in different nights under 
varying photometric and seeing conditions. Standard star 
observations revealed that the conditions were not pho- 
tometric during all nights and an adjustment (described 
below) was necessary to align all nights onto a photomet- 
ric system. As in Papers I and II we will always refer to 
the Kg filter as K. 

The NGC 5846 data were obtained during the nights 
of April 6th, 8th, and 9th 1999 while the NGC 7192 
data were taken during the nights of June 1st, 2nd, 21st, 
and 22nd. The observing strategy for NGC 5846 data 
was the following: lOxlOsec object -|- 5x(2xl0sec) sky -I- 



Hempel et al.: Extragalactic Globular Clusters in the Near Infrared III 



3 




Fig. 1. Field of view of NGC 5846 (upper) and NGC 7192 (lower). The ISAAC field was chosen to fit the available 
archive data taken with WFPC2 onboard the HST. The dashed lines in both frames indicate the HST field and the 
solid one the corresponding ISAAC image. Dotted circles show the 4 and 8 kpc distance to the galactic centers, the solid 
circles the half light radii of the galaxy's light. The globular cluster samples were divided in two colour sub-populations. 
Dots mark objects with 0.8<{V — I)<1.2 and open circles represent clusters in the range 1.2<(I^ — /)<1.5. 



lOxlOsec object. The one for NGC 7192: lOxlOsec object 
+ 2x(6xl0sec) sky + lOxlOsec object. The monitoring of 
the seeing, magnitude variations, sky level and geometric 
correction was done by measuring the parameters of four 
isolated stars in each single frame for both galaxies. A 
detailed description of the sky subtraction and combina- 
tion procedure of the K band data is given in Puzia et al. 
(2002). The FWHM of the stellar PSF in the final iC-band 
image is -0.4" for NGC 5846 and -0.5" for NGC 7192. 
The total exposure times for NGC 5846 and NGC 7192 
are 10000 sec, and 12000 sec, respectively. 

2.1.1. Photometry 

The photometric calibration of the NGC 5846 and 
NGC 7192 data set was based on the photome- 
try of 3 and 4 near-IR standard stars, respectively 
UPersson et al. 1998|l . 

Measuring the instrumental magnitudes in an aperture 
of 6 pixel diameter and applying the same analysis as de- 
scribed in Puzia et al. (2002), the following calibration 
relations for the photometric nights have been derived 

i^5846 = hnst + 23.81(±0.018) - 0.05(±0.009)x (1) 
K7192 = fcinst + 23.83(±0.010) - 0.05(±0.007)x (2) 

where -ft'gaiaxy is the calibrated magnitude, fcjnst is the in- 
strumental magnitude, and x the effective airmass {— 1.23 
for NGC 5846 and 1.34 for NGC 7192). The error of the 
zero points (second term in eq.(l) and (2)) includes pho- 
tometric errors of each single standard star measurement 
and the errors of the aperture correction analysis. The er- 



ror of the airmass term is an estimate from the variations 
in airmass of all single exposures. 

The zero point shifts of the non-photometric nights to 
photometric conditions has been derived by tracing the 
magnitudes of four isolated stars over all nights. The final 
photometry was performed on the overall combined image. 
For NGC 7192 it was first done on two combined images 
separately. One for the photometric nights (June 1st and 
2nd), and one for the non-photometric nights (June 21st 
and 22nd) . The magnitudes of both images have been av- 
eraged after a zero point correction of the latter. The true 
photometric uncertainty, measured by the scatter of the 
single measurements is of the order of <0.04 mag, mainly 
due to the strongly varying sky background. 

Finally, all magnitudes were corrected for Galactic 
foreground reddening using the reddening values of Table 
1 and the extinction curves of Cardelli et al. (1989). The 
corrections for NGC 5846 and NGC 7192 are Ak = 0.020 
mag and Ak = 0.012 mag. 

2.2. HST/WFPC2 Optical Data 

The HST data were taken from the public HST archive. 
NGC 5846 was observed with HST -I- WFPC2 under pro- 
gram GO.5920. The total exposure times of the combined 
images are 6600 sec in F555W, and 6900 sec in F814W. 
NGC 7192 has been imaged with WFPC2 under program 
GO. 5943 in F555W and F814W filters with 1300 sec and 
1000 sec of total exposure time, respectively. The HST 
images were reduced and calibrated following the proce- 
dure as described in Puzia et al. (1999, 2002). All mag- 
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nitudes were measured with the SExtractor tool using 
a 8 and 4-pixel-diameter aperture for the PC and WF 
chip and corrected with respect to the Holtzman 0.5" 
standard aperture HHoitzman et al. 1995|> . Instrumental 
magnitudes were then transformed to the Johnson V 
and / magnitudes according to the prescription given by 
Holtzman et al. (1995). All magnitudes were reddening 
corrected using the following values: Ay = 0.182 and 
Ai = 0.107 for NGC 5846 and Ay = 0.113 and A/ = 0.066 
for NGC 7192 (see Eb-v in Table 1). 

2.3. Selection criteria 

After combining the optical and near-infrared data us- 
ing the GEOMAP task within IRAF, the globular cluster 
sample includes 184 and 61 objects for NGC 5846 and 
NGC 7192, respectively. Figure H shows the field of view 
for both galaxies and the spatial distribution of the cluster 
candidates. In order to limit the contamination of the sam- 
ple by background galaxies or foreground stars and to set 
a limit onto the photometric error, general selection crite- 
ria have been applied. In our discussion only objects with 
an photometric error A(y — /) < 0.1 mag and A{V — K) < 
0.1 mag and a FWHM of the PSF below 0.25" in the V 
and I-band are considered. 

3. Colour-colour diagrams for NGC 5846 and 
NGC 7192 

Colour-colour diagrams together with various SSP mod- 
els (e.g. in this paper by Bruzual & Chariot 2000) are the 
basis for age and metallicity estimates. Different models, 
at identical colors, can show differences in absolute age of 
about 3 Gyr. Thus, the SSP approach can only lead to ap- 
proximate absolute ages for the sub-populations. However, 
relative ages can be estimated accurately enough to sepa- 
rate sub-populations built up during major star formation 
events. 

As shown in Puzia et al.(2002), the various models dif- 
fer mainly in the metal-rich range (approx. [Fe/H]> —0.8). 
A detailed comparison between the different models in an 
absolute sense can be found in Maraston et al. (2001a). 
The metal-rich regime, however, is exactly the one we in- 
tend to probe, given that the intermediate-age populations 
are expected to be enriched in metals. 

As opposed to the situation for absolute ages, the rel- 
ative age predictions for given colors are relatively similar 
from model to model. Since we focus on the relative age 
dating of globular cluster systems the choice of a specific 
model is not crucial and we use the model of Bruzual & 
Chariot (2000) throughout the following analysis. 

As we will show in the following two subsections the 
results of the KMM test l|McLachlan fc BasfordlMsl up- 
dated version 2001) and the derived metallicities are found 
to differ significantly from expected values and even seem 
to be in contradiction to what visual inspection of the 
colour-colour diagrams tells us. Considering the relatively 
small sample size, the photometric errors in (V — K) and 



the limited depth of the observations (shifting the peak 
position toward red colors) the results of the KMM test, 
as described by McLachlan & Basford (1988), have to be 
discussed with much caution. Since this work concentrates 
on relative ages at this stage we will only mention the re- 
sults of the KMM and leave the discussion for later. 

3.1. NGC 5846 

The colour-magnitude diagram (hereafter CMD) and 
the colour-colour diagrams for NGC 5846 are given in 
Figures El and El The histogram in the upper part of the 
CMD shows the complete set of clusters as a solid his- 
togram whereas the selected objects, following the selec- 
tion criteria given in Section lT^ are shown by the hatched 
histogram. As expected the selection criteria affect mostly 
the objects in the very red colour range, since the photo- 
metric errors are larger and background galaxies possibly 
contaminate the sample. 

The CMD (Fig. O reveals evidence for the a bi- 
modal colour distribution. The KMM test as described 
by Ashman et al. (1994) confirms this result with a con- 
fidence level of ~90% for a bimodal colour distribution. 
We obtain peak positions of {V — -ft')biuc = 2.57±0.06 and 
{V — K)i-cd — 3.18±0.06 for both colour populations with 
an estimate of correct allocation value (confidence level) of 
0.79 and 0.94 for the blue and the red peak, respectively. 
Hereby about 30% of the globular clusters were assigned 
to the blue population and 70% to the red. Using the 
calibration by Kissler-Patig (Paper I) for old populations, 
this corresponds to a metallicity of [Fe/H]=— 0.54±0.6 dex 
and 0.27 ± 0.4 dex for the metal-poor and the mctal-rich 
populations, respectively. 

Note that our completeness limit is dominated by the 
i^-band, and that we are biased in favour of red clusters 
(Paper II; see Fig.j^J. Thus, both the ratio of the numbers 
of red to blue clusters and the metallicity of the blue peak 
are overestimates. 

It is interesting to mention that no multi-modality 
could be found in the V — I colour distribution. This 
is similar to the situation in NGC 4365 where the 
intermediate-age population, when projected on the V — I 
axis, fills the gap between the two old (metal-rich and 
metal-poor) sub-populations. When forced to bimodality, 
the {V — I) peak positions determined by the KMM test 
are {V ~ /)biuo = 1-11 and {V - I)rcd = 1-13 ± 0.02. 
This is in agreement with the results by Gebhardt & 
Kissler-Patig (1999). 



3.2. NGC 7192 

Referring to the CMD and the colour-colour diagram for 
NGC 7192 given in Figure^and FigurejSl the difference to 
the globular cluster system of NGC 5846 can be seen: while 
NGC 5846 exhibits a strong spread perpendicular to the 
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(V K)„ [mag] 

Fig. 2. V vs. {V — K) colour-magnitude diagram for NGC 5846. The top sub-panel shows the colour distribution of all 
(open histogram) and selected (hashed histogram) objects. The solid line marks the probability density distribution 
together with its la uncertainty (dotted line). The lower sub-panel shows the CMD. Here the filled symbols mark 
the selected clusters while open circles indicate rejected objects (see Section |^3J • Since the photometric errors are 
dominated by the X-band, only the {V — K) errors are shown. The dashed line marks the limiting magnitude in the 
K-himd {K = 21.5 mag). 
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isochrones, NGC 7192 is more homogeneously populated 
along the isochrone. 

The small number of objects do not allow one to draw 
strong conclusions from the KMM test in terms of the 
existence of colour bimodality. The colour-colour plot for 
NGC 7192 does suggest the existence of two populations of 
different metallicity (i.e., a 'metal-poor' population with 
{V -I) < 1.1 and {V -K) < 2.8 and a 'metal-rich' pop- 
ulation with {V - I) > 1.0 and {V - K) > 2.8) but the 
small number statistics do not allow one to support this 
firmly statistically. 

Formally, the colour for the two peaks of the distribu- 
tion projected on the V — K axis is returned by KMM 
to be {V - K)= 2.83 and 2.89 (with 90% of the objects 
assigned to the first peak), if a bimodal distribution is 
assumed. Since the difference in colour between the two 
peaks is well below the photometric error it seems far- 
fetched to assume the existence of two distinct globular 
cluster sub-populations. 

Therefore we assume that the NGC 7192 globular 
cluster system consists only of one dominant population. 
Following the calibration values by Kissler-Patig (Paper 
I) we derive a peak metallicity of [Fe/H]=— 0.57 ± 0.37. 

4. Determining the age structure from 
colour-colour diagrams 

4.1. Colour distributions 

The basic idea in using optical and ncar-infrarcd photom- 
etry in globular cluster studies is to lift the age - metallic- 
ity degeneracy and to resolve possible age sub-populations 
of globular clusters. The minor drawback with respect to 
using optical colors alone is an increased photometric er- 
ror mainly caused by the infrared observations, i.e., in the 
K-band. This, however, is more than compensated for by 
the larger diagnostic power of the optical - near-infrared 
combination (see Papers I and II). A separation in the 
colour-colour plot of sub-populations with different metal- 
licities becomes easier, and a separation in age at a given 
metallicity becomes feasible. 

We investigate below what the diagnostic power of 
the method is in terms of separating sub-populations of 
different ages. To do so, we investigated artificial colour- 
colour distributions of composite populations based on the 
Bruzual & Chariot (2000) isochrones. These artificial sys- 
tems were built as follows: 

• In a first step we create an artificial V — K distribu- 
tion of metal-rich globular clusters. We populate the red 
{V—K) colour range (2.7 < {V—K) < 3.8), assuming that 
it is occupied by "old" and "intermediate" age objects. 
We do not consider bluer clusters (2.0 < (V - K) < 2.7), 
assuming that these are only "old" objects. 

For now, we assume cases similar to NGC 5846 (see 
Figure ISJ in terms of sample size and photometric er- 
rors. For a first exemplary case, we assume 50% of the 
red population to be 15 Gyr old and 50% to be 3 Gyr 
old. These numbers are not assumed to reflect the 'true' 



situation, but rather serve to demonstrate the method at 
this point. We will probe models with different ratios in 
the future (see Section IBT^ . Further, it will become clear 
below that we are not directly comparing observed and 
simulated colour-colour diagrams but rather their cumu- 
lative age distributions, i.e., properties of the distributions 
still need to be "calibrated" . 

The final modeled systems contain 43 old, blue ob- 
jects (not considered further) and 120 red objects homo- 
geneously distributed within the V — K range ^ . The red 
population was divided into an old (15 Gyr) and young (3 
Gyr) population with 60 objects each. 

• In a second step we attach to each V — K data 
point a (Icr) error drawn randomly from our observed list 
of A{V — K) errors for NGC 5846, and then smear in a 
Monte-Carlo approach each V — K point with up to ±3 
times its associated error (i.e. allowing up to 3a errors in 
very rare cases) . These new V — K values are stored with 
their 1 a error and used for the further process. 

• The third step consists of associating a.V — 1 colour 
to each new V — K colour. To do this, we perform a least- 
square fit to the SSP model isochrones (in this case from 
Bruzual & Chariot 2000) by a logarithmic function {{V — 
I)=A-\og{V — K)+B) . The particular fit parameters A and 
B are given in Table 2. Figure shows the isochrones as 
given byp3ruzual 2000 and our fits to the isochrones (solid 
lines). 

The fits are then used to compute for each V — K 
point the corresponding V — I point, once the age was 
chosen. In our case, we used the 15 Gyr fit for the 60 
old artificial clusters, and the 3 Gyr fit for the other 60 
(young) artificial clusters. 

• Finally, in the fourth step we associate a mea- 
surement error A(V^ — /) to each {V — I) data point in 
a similar way as for {V — K). We now have a set of 120 
objects (60 old, 60 young) with associated {V — K) and 
{V — I) colors and errors. 

As an example of such artificial colour-colour distri- 
butions, we show in Figure [71 the modeled colour-colour 
diagram of a purely 15 Gyr old population and a com- 
posite old (15 Gyr) and young (3 Gyr) clusters as well 
as the resulting age distribution (lower panels). This dis- 
tributions should be compared to the observed data for 
NGC 5846 in Fig. El 

4.2. Cumulative age structure 

The separation of the isochrones in Fig. and the differ- 
ence between the distributions in Fig. [7\ (upper panels) 
motivated us to use the optical-near infrared photometry 
in order to try to separate globular cluster populations of 

^ We experimented also with Gaussian distributions within 
this colour range, but the final results of the experiment did not 
differ significantly and we have no better physical justification 
for a gaussian than for a homogeneous colour distribution 
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Fig. 4. V vs. iy ~ K) colour magnitude diagram for NGC 7192. (panels and symbols as in Fig. [2J. 
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Fig. 5. Colour-colour diagram for NGC 7192. No clear indication for a second red sub-population can be found 
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Table 2. Fit parameter for model isochrones IjBruzual 2000|l for a fogarithmic fit {V - /)= A-fog(y - K)+B. 



Isochrone [Gyr] 


f 


2 


3 


5 


7 


10 


13 


15 


A 


0.5645 


0.7136 


0.7674 


0.7795 


0.7922 


0.8134 


0.8423 


0.8573 


B 


0.3296 


0.2809 


0.2633 


0.2765 


0.2826 


0.2766 


0.2596 


0.2502 




Fig. 6. SSP isochrone fit for {V - /) vs. {V - K) colom- 
colour diagrams. The symbols are colours given by the 
SSP models fBr uzual 2000|l for 1 Gyr (open triangle), 3 
Gyr (sohd triangle), 5 Gyr (open square), 10 Gyr (solid 
square) and 15 Gyr (dots). The lines represent the result of 
a logarithmic fit to the isochrones with {V—I)= A-log(T^— 
K) + B. The fit parameters A and B are given in Table 2. 



Fig. 7. The upper two panels show the {V — I) vs. {V — K) 
colour-colour diagram for a single age population (left 
panel) and a composite 50% old (15 Gyr) and 50% young 
(3 Gyr) cluster population (right panel). In the lower pan- 
els, the cumulative age structure (see ^4.2|l for both cases 
is given as the mean of 1000 such simulations. The 50%- 
level is marked as a solid line. 



different ages. This appears from the above plots to be at 
least possible to age differences of 10 Gyr or more (with 
respect to 15 Gyr for the oldest clusters). Therefore we 
decided to take the analysis one step further and to make 
a first attempt towards determining the age structure of 
the globular cluster systems. 

We first attempted to compare directly the observed 
and the artificial colour-colour distributions, but this 
turned out to be too complicated. Two dimensional statis- 
tical tests are required in that case and the typical number 
of data points is not leading to any statistically signifi- 
cant results. We experimented with 2-dimensional f-test, 
2-dimensional Kolmogorov-Smirnov tests, but in order to 
constrain the tests somewhat, one needs to put rather ar- 
tificial constraints and the physical meaning of the final 
results is dubious. 

We therefore decided to focus on 1-dimcnsional repre- 
sentations of the age structure. The dimension was chosen 
to maximize the age gradient along it, i.e., "perpendicu- 
lar" to the isochrones. 



Briefly, we associate to every cluster in our observed 
or artificial distribution an "age greater than X" when 
it lies above the isochrone of that age X in the colour- 
colour diagram. We start with the youngest isochrone 
(0 or 3 Gyr, see below) for which most cluster will lie 
above, and then move up isochrone by isochrone (in the 
steps 1,2,3,5,7,10,13,15 Gyr). Hereby the notation "0 Gyr" 
refers to objects below the 1 Gyr isochrone. The result is 
an inverted cumulative distributions as shown in Fig. [S] 
for artificial distributions and in Fig. for the observed 
distributions of all our galaxies analyzed so far. The cu- 
mulative age distribution can be represented in absolute 
numbers (left panel) or normalized to the total number of 
objects (at an arbitrary bin, right panel). 

4.2.1. Artificial data sets 

As described above, the first set of simulations was done 
for combinations of a 15 Gyr old sub-population and 
an equal number of intermediate-age objects (1,2,3,5,7,10 
Gyr). The results are shown in Figure|Hl Each curve shows 
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Fig. 8. Age distribution for mixed cluster populations as- 
suming half of the objects being 15 Gyr old and half of 
them being respectively 1, 2, 3, 5, 7 or 10 Gyr old. The 
50% level is marked by a solid line. 



the mean age structure as evaluated from 1000 models for 
that given age composition. 

The realistic photometric errors create a spread in the 
colour-colour diagram such that even a pure 15 Gyr sys- 
tem does not show 100% of the clusters to be older than 
15 Gyr. Instead, the spread in the {V — I) vs. {V — K) 
diagram leads to a gentle fall-off with isochrone age. 

However, this fall-off clearly changes when a second, 
younger sub-population gets mixed in. By the time one 
mixes a 1,2,3 Gyr population in the system, the fall-off 
becomes very steep around 2,3,4 Gyr, and the curves cross 
the 50% level well before the 10 Gyr mark. 

4.2.2. Observed data sets 

The result for the artificial distributions can be compared 
with the observed age distributions of the globular clus- 
ter systems (see Figure in galaxies analyzed so far 
(M 87, NGC 4478, NGC 4365, NGC 3115, NGC 5846 and 
NGC 7192). 

Taking first the observed data alone, we notice a clear 
similarity between NGC 5846 and NGC 4365, as opposed 
to the 4 other systems (see right panel of Fig. |51 Both 
galaxies fall-off steeply at early isochrone age and cross 
the 50% line well before 10 Gyr. 

When compared to the results of artificial distribu- 
tions, this leads immediately to the interpretation that 
both must host a significant fraction of intermediate-age 
clusters within their red sub-population. For NGC 4365, 
this was suspected already in Paper II and has since then 
been confirmed spectroscopically (Larsen et al. 2002). 



In contrast, the age distribution of NGC 7192, 
NGC 3115, NGC 4478 and M 87 seems to be more consis- 
tent with what we would expect for a single age and old 
(>10 Gyr) population. 

4.3. Contamination of background objects 

Contamination of background objects is a potential prob- 
lem in our analysis and we briefly investigate its impact 
below. 

As above for the artificial distributions, we adopt 
a situation similar to the observations of NGC 5846. 
We used the Hubble Deep Field South (available at: 
www.stecf.org/hstprogrammes/ISAAC), which presents 
two advantages: it covers exactly the same (WFPC2) field 
of view as our observation, and it has a deep enough K- 
band observation to match our ISAAC observations of the 
globular cluster systems {K < 21.5 mag). The HDF-S 
sample was further selected in colours as for our samples. 
However, we could not apply a FWHM selection (having 
only the list of objects) so that the contamination is ex- 
pected to be an overestimate. Even in that case, however, 
we show below that the effect is negligible for cases such 
as NGC 5846. 

FigureEHshows the colour-colour diagram of the HDF- 
S. Depending on magnitude selection, we have between 
25 and 40 background objects in our colour selection box. 
The resulting cumulative age distribution for the HDF-S 
sample is shown in Figure ^2 and assigns about 50% of 
the objects to a population younger than 3 Gyr. Further, 
from the colour-colour plots, it becomes clear that the 
majority of these objects actually lie below the 2 Gyr line, 
being bluer in {V — K) and {V — I) than the intermediate- 
age, metal-rich globular clusters in NGC 5846. It is most 
probably dominated by a star-forming galaxy population 
that would be rejected by our FWHM criteria. 

The effect of contamination on NGC 5846 is shown 
in Figure El There, we plot the uncorrected (absolute 
and relative) age distributions, as well as the ones cor- 
rected for background contamination using the HDF-S 
source counts. Clearly, the effect is marginal. The con- 
taminating objects tend to drag the age distribution to- 
wards younger ages, mimicking a slightly younger sub- 
population, e.g. the intersection with the 50%-level oc- 
curs at a larger age. Thus, "old" distribution such as for 
M 87, NGC 3115 etc would appear even older when cor- 
rected for contamination. In the cases of NGC 4365 and 
NGC 5846 (large numbers of clusters) the effect of this 
(overestimated) contamination is small and does not in- 
fluence the conclusion that these systems host a significant 
intermediate-age population. 

5. Conclusions and future work 

5.1. Conclusions 

In this paper, we presented results of our optical - near- 
infrared study of the globular cluster systems of NGC 5846 
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Fig. 9. Age distribution for different galaxies (see text). The absolute number counts (left panel) are normalised to 
the total numbers of clusters in the sample (right panel). It is clearly seen that two systems (NGC 5846, NGC 4365) 
are significantly different from NGC 4478, M 87, NGC 7192 and NGC 3115. This is interpreted as NGC 5846 and 
NGC 4365 hosting intermediate-age sub-population (which was confirmed spectroscopically for NGC 4365). The 50% 
level is marked by a solid line. 




2.5 3 

(V-K)„ [mag] 



2.5 3 

(V-K)„ [mag] 



Fig. 10. {V — I) vs. {V — K) colour-colour diagram for the 
Hubble Deep Field-South using different K completeness 
limits: left, all objects with K < 21.5, right all objects with 
K < 21.0 (bracketing our completeness limit for NC 5846). 
The box marks the colour range in {V — I) and {V — 
K) used for the determination of the globular cluster age 
structure as described in section Both diagrams show 
that the highest contamination of our sample is expected 
below the 2 Gyr isochrone. The isochrones superimposed 
are from Bruzual & Chariot SSP models (2000). 

and NGC 7192. While the latter does not show any sig- 
nificant anomaly, the former shows a case very similar to 
that of NGC 4365 (see Paper II) with good evidence for a 
significant metal-rich, intermediate-age population in ad- 
dition to an old, metal-rich one. 

We developed a new method to quantify the age struc- 
ture of the globular cluster systems studied so far, based 
on a comparison of the cumulative age distributions of 
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Fig. 11. Age structure of the HDF-South objects using the 
procedure given in Section f4.2l Left and right panels show 
the absolute and relative age distributions respectively. 
The colour selected samples with K < 21.5 and K < 21.0 
are shown as solid and dashed curves, respectively. The 
50% level is marked by a solid line. 



the observed and modeled optical - near-infrared colour- 
colour diagrams. This method appears powerful enough 
to detect intermediate-age sub-population within globu- 
lar cluster systems. 

Our first conclusions for the galaxies studied so far are 
that both NGC 4365 and NGC 5846 have a cumulative 
age distribution of the metal-rich globular clusters that is 
different from those of the 4 other galaxies (NGC 3115, 
NGC 4478, NGC 7192 and JVI 87). The former are better 
modeled by by a composite metal-rich population includ- 
ing an old (15 Gyr) as well as a significant young (1-5 
Gyr) population. The others, in contrast, are very similar 
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Fig. 12. Influence of background contamination on the 
age structure of the system in NGC 5846. The left and 
right panels show the absolute and relative cumulative age 
distributions, respectively. The uncorrected distribution 
is shown as a solid curve. The distribution corrected for 
background contamination of K < 21.5 and K < 21.0 are 
shown as long dashed and short dashed lines, respectively. 

to the model for pure old populations and are thus best 
explained by being dominated by old objects. 

NGC 7192 suffers from small number statistics, but 
the current dataset is, with respect to the colour-colour 
diagram and the age structure, more consistent with only 
one sub-population. Compared to the globular cluster sys- 
tems investigated so far, NGC 7192 most closely resembles 
NGC 3115. 

For M 87 and NGC 4478 the results given in Paper 
I were not fully conclusive. NGC 4478 also suffers from 
small number statistics. M 87 appears now more clearly 
dominated by two sub-populations differing in metallic- 
ity, but with a metal-rich population dominated by old 
objects (see also Jordan et al. 2002). 

We need to emphasize again that the method cannot, 
currently, produce reliable absolute numbers, neither in 
terms of ages, nor in terms of ratios between the different 
populations. The results in this paper should therefore 
be considered as qualitative for now, and will be better 
quantified in the future. 

5.2. Future work on age dating 

The upcoming papers will analyze the remaining galaxies 
in our dataset (in total 11 galaxies). This will allow us 
to discuss the results in the light of galaxy properties. In 
particular, we will look at trends with galaxy size and 
environment. The local density (TuUy 1988) ) covered by 
our galaxy sample spans from 0.08 Mpc~'^ (NGC 3115) 
to 4.17 Mpc~^ (M87) and environmental effects appear 
to be an important ingredient to galaxy formation and 
evolution. 

On the modeling side, our goal is to improve the quan- 
titative information of our method. We are developing a 

test to find the best solution in the two parameter space 
of the models: ratio old/young and age of the young pop- 



ulation, which are slightly degenerate. We also plan, with 
the help of spectroscopy and wide field photometry, to be 
able to better calibrate the models in terms of absolute 
age. 

Also, we are exploring the dependence of the results on 
a particular SSP model and we will repeat the determina- 
tion of the age structure using SSP models by Maraston 
(2000) and Vazdekis (1999). Further, we will investigate in 
more detail the effect of background contamination. This 
latter aspects will be the subject of a separate paper. 
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